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Abstract 

We employ the linear cr model to study the chiral dynamics of two flavors of quarks at finite 
temperature and density. Our calculations include thermal fluctuations of both the bosonic and 
fermionic fields. In particular, we determine the phase diagram in the plane of temperature and 
baryon chemical potential as a function of the pion mass. An interesting phase structure occurs 
that results in zero, one, or two critical points depending on the value of the vacuum pion mass. 



The physical pion mass is small but not zero. In consequence, the conventional wisdom is that 
there is no true thermodynamic chiral phase transition at finite temperature T and zero baryon 
chemical potential yU. Instead, there is expected to be a curve of first-order phase transition in 
the pt-T plane that terminates in a second-order phase transition at some critical point (jj. c , T c ). 
The location of the critical point obviously depends on the physical (vacuum) pion mass. This 
topic has been under intense theoretical study using various effective field theory models, such 
as the Namu Jona-Lasinio model JH S 0], a composite operator model J3l, a random matrix 
model [5], a linear cr model J3l, an effective potential mod el Iftll, and a hadronic bootstrap model 
[7], as well as various implementations of lattice QCD Iill^ ll0l[l"lll . It is also of great interest 
because collisions between heavy nuclei at medium to high energy, such as at the future Facility 
for Antiproton and Ion Research (FAIR), or possible low energy runs at the Relativistic Heavy 
Ion Collider (RHIC), may provide experimental information on the phase diagram in the vicinity 
of a critical point. In this work we study the phase diagram of the linear cr model coupled to 
two flavors of identical mass quarks. Whereas the cr model is an oft-used effective model that 
represents some of the essential features of the chiral dynamics of QCD, the reason to couple 
the fields to quarks is less obvious. One argument is based on the existence of the critical point 
itself. If a critical point exists, then one can go around it without crossing the curve of the 
first-order phase transition. The effective degrees of freedom should not change too much in 
following such a path. Therefore, if quarks are considered to be reasonably useful degrees of 
freedom on the higher temperature side then they should be useful on the lower temperature side 
too. For this reason we use quarks to carry the baryon number but acknowledge the resulting 
uncertainty in the results. The present work extends that of Refs. O O, H Q in several 
ways. We include thermal fluctuations of the meson and fermion fields, which can be important 
at finite temperature when the magnitude of the fluctuations becomes comparable to or greater 
than the mean values. We scan vacuum pion masses from zero to over 300 MeV, which may be 
particularly useful for comparison with lattice gauge theory with different quark masses. Even 
when all other parameters of the model are fixed, we find an interesting phase diagram that may 
have zero, one or two critical points depending on the value of the vacuum pion mass. 
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The Lagrangian is 



1 9 1 2 

£=-(d l j,7r) +-(<Vr) - U(o-,7i) + i(r[ifl-g(o- + iy 5 T-7i)]i[/, (1) 

where 

U(cr, 7t) = ^ (a 2 + n 2 - f 2 f - Her (2) 

in an obvious notation. The SU(2)iXSU(2)« chiral symmetry is explicitly broken by the term 
Her, which gives the pion a mass. The parameters in the Lagrangian are constrained by fixing 
the pion decay constant f n = 92.4 MeV, the er mass m a = 700 MeV, and the quark mass is set 
to one-third of the nucleon mass m q - 313 MeV. The vacuum pion mass m„ is varied from to 
nia-12 = 350 MeV. To create an effective mesonic model, we integrate out the quark degrees of 
freedom in the usual way flill . such that 

InZquark = IndetD, (3) 

where D is the inverse quark propagator. This leads to an effective Lagrangian 

Xeff = \ (d^f + l - (d^erf - U eB (cr, Jt) , (4) 

where 

T 

U eS {cr, 7T) = U (cr, n) - — In Zquark(cr, 7r) (5) 

is the effective potential. We decompose the scalar field into a condensate v plus a fluctuation A. 
The thermodynamic potential is computed from 

A = <^eff> - \ml(A 2 ) - l -m 2 n {7! 2 ) + Qtr + £2„ , (6) 

where Q., T and Q.^ are the independent particle contributions from the cr and pion quasiparti- 
cles. Angular brackets refer to thermal averaging. The condensate and quasiparticle masses are 
determined self-consistently from 

du e A 2 ld 2 u e A 2 ld 2 u e A 

«V = ttt" = . (7) 



dv I " \ 8A 2 I * \ dn 

The resulting equations must be solved (very accurately) using numerical methods. All thermo- 
dynamic identities were verified both analytically and numerically. 

An example of the pressure versus chemical potential at fixed temperature is shown in the left 
panel of Fig. 1. For the higher temperature of 80 MeV there is only one self-consistent solution. 
For the lower temperature of 50 MeV, there is a unique solution at large yU and another unique 
solution at small pi. For a range of yU centered about 900 MeV there are three solutions: One is 
associated with a continuation of the low-density phase, one is associated with a continuation 
with the high-density phase, and the third solution (not shown in the figure) is an unstable phase. 
The point where the two curves cross is the location of the phase transition. In this example it 
is first-order since there is a discontinuity in the slope dP(jj,T)/dfi = hb- Where each curve 
terminates is the limit of metastability for that phase. The thermodynamically favored phase 
is the one with the largest pressure. At some temperature between 50 and 80 MeV the slopes 
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Figure 1: Left panel: Pressure versus baryon chemical potential for temperatures below and above T c for the physical 
vacuum pion mass. Right panel: Phase diagram for illustrative values of the vacuum pion mass. 



are equal at the crossing point, there are no metastable phases, and the second derivative is 
discontinuous. This is indicative of a second-order phase transition. The location of this point in 
the fi-T plane is the critical point. 

The phase diagram for a sampling of vacuum pion masses is shown in the right panel of Fig. 
1 . Generically there is a curve of first-order phase transition starting on the // axis and arching 
to the left. This curve terminates at a critical point. For m„ VSLC = 321 MeV the critical point sits 
on the T — axis, and for m^ vac > 321 MeV there is no phase transition at all. Of course the 
precise numbers depend on the constants in this model, such as the vacuum cr mass and what 
value one assigns to the coupling of the quark field to the <x field, but the results are in line with 
expectations II 1 711 . 

Something very interesting happens as the vacuum pion mass is decreased from 50 MeV to 
0. The left panel of Fig. 2 shows the phase diagram for a vacuum pion mass of 35 MeV. There 
are now two critical points! There is a line of first-order phase transition beginning on the n axis 
and arching to the left to end at a critical point of fi c i * 725 MeV and T c \ » 92 MeV and there 
is another line of first-order phase transition beginning on the T axis and arching to the right to 
end at a critical point of // C 2 ~ 240 MeV and T C 2 ~ 137 MeV. The latent heat curves are shown 
in the right panel of Fig. 2. In the vicinity of fi — 575 MeV and T — 1 10 MeV the latent heat 
gets pinched to zero for some critical value of the vacuum pion mass between and 35 MeV. 

The results shown in the right panel of Fig. 1 illustrate the conventional view on the nature 
and location of the critical point. The result shown in the left panel of Fig. 2 is unconventional 
or exotic. To the best of our knowledge, this is the first time a model calculation has resulted in 
such a phase diagram. It cannot be expected that a model as simple as the one studied here can 
make definitive predictions for what actually happens in QCD with the observed values of the 
quark masses. However, it can serve to illustrate the possibilities. 

More details on this work have been, or will be, presented elsewhere |18, 3. Extensions of 
it can be easily identified and are under current investigation. 
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Figure 2: Left panel: Phase diagram for a vacuum pion mass of 35 MeV. There are two critical points. Right panel: 
Latent heat along the coexistence curve for illustrative values of the vacuum pion mass. Note the second region in the 
lower right corner for the 35 MeV mass. 
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